Uniform curves for van der Waals interaction between single-wall carbon nanotubes 
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We report very simple and accurate algebraic expressions for the van der Waals (vdW) potentials 
and the forces between two parallel and crossed carbon nanotubes. The Lennard- Jones potential for 
two carbon atoms and the method of the smeared out approximation suggested by L.A. Girifalco 
were used. It is found that interaction between parallel and crossed tubes are described by different 
uniform curves which depend only on dimensionless distance. The explicit functions for equilibrium 
vdW distances, well depths, and maximal attractive forces have been given. These results may be 
used as a guide for analysis of experimental data to investigate interaction between nanotubes of 
various natures. 
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The van der Waals (vdW) interaction between 
graphitic structures is very important for application in 
Nano Electro Mechanical Systems. There are a number 
of publications devoted to estimation of vdW potentials 
for graphite layersi^S.^ two fuUerenes^i^i^i^, fuUerene and 
surfaceii^, carbon nanotube (CNT) and surface^ i^°'^^ , 
fuUerenes inside and outside of nanotube s ^^i-^^'-^^i-^^'-^^ . 
The interactions between the inner and the outer paral- 
lel tubes such as single- (SWNTs>iaii2a8ii9, double-^OiSi^ 
and multi-wall nanotubes (MWNTs) ^^i^^'^'* are also well 
studied. The potential between two crossed CNTs was 
discussed in our previous work^^. 

The continuum Lenard- Jones (LJ) model suggested by 
L.A. Girifalco^ is usually used to evaluate potential be- 
tween two graphitic structures. The LJ potential for two 
carbon atoms in graphene-graphene structure is 
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where r is a distance, A and B are the attractive and 
repulsive constants. 

The potential between two SWNTs is approximated 
by integration of LJ potential 
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where dSi and are the surface elements for each 
tube. In the case of CNT-CNT interaction the mean 
surface density of carbon atoms is = 4/V3a^, where 
a = 2.46 [A] is the lattice constant for graphene hexago- 
nal structure. 

The LJ potential from Eq. (2) can be integrated ex- 
actly for two crossed CNTs^. Unfortunately the ana- 
lytical formula is very complicated expression in terms 
of elementary functions and elliptic integrals. In the 
case of parallel tubes the numerical integrations were 
appliedii^'SS. 

It was found that the vdW potential ip between Cgo — 
Ceo, Cgo-SWNT, Cgo-graphene, graphene-graphene, and 
parallel SWNT-SWNT or MWNT-MWNT at different 
distance d can be described by the universal curve^^ d'^'i^'^ . 



The universal curve for two tubes means that a plot of 
Tp = ip/\ipQ\ against d = djd^ gives the same curve for all 
radii of tubes, where (^o is the minimum energy and dg is 
the equilibrium spacing for the two interacting surfaces. 

In the present work we report very simple and accurate 
algebraic formulas for vdW potentials and forces between 
parallel and crossed SWNTs. We declare that the inter- 
action between parallel and crossed tubes are described 
by different uniform curves. Also we give explicit func- 
tions for equilibrium vdW distances, potential wells, and 
maximal attractive forces. 

Figure 1 illustrates the interaction between two 
SWNTs. In this figure t\ and t2 are the radii, R is the 
distance between axes, and d is the gap between surfaces 
of tubes. If angle 7 = then tubes are in parallel. 

Let's consider first the case of crossed tubes. In our 
previous work^- we obtained for two crossed SWNTs that 
the vdW potential is 
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where gA and gs are complicated expressions in terms of 
elliptic functions. 

These expressions can be essentially simplified when 
d/ti <C 1 and d/t2 ^ 1. Introducing a = 1/ti + 1/^2 
and using expansion in small parameter we obtain the 
approximating formulas 
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This approximation allows us to explain the existence 
of uniform curves for crossed SWNTs. Also it allows to 
get simple expressions of equilibrium vdW distance do, 
minimum vdW potential (^Qj uniform curve ^ = (/3/|iy9o|i 
and total force F for interaction between two crossed 
SWNTs. 

From Eq. (3) using Eqs. (4) and (5) we can find 
dip/dd = and write the recurrent equation for equi- 
librium distance cIq: 
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One can solve Eq. (6), passing a few steps of this 
recursive formula with very fast convergence. 

If ado 1, that is ti and t2 are big enough in com- 
parison with do, then we have the first approximation for 
equilibrium spacing 
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Remarkable that for SWNTs of large radii, the equi- 
librium distance depends only from the attractive and 
repulsive constants A and B. Measuring the gap do it is 
possible to find the ratio between A and B. 

We find that Eq. (7) works weU if ii,i2 > 6 [A] and 
gives small error for tubes of smallest radii. The second 
approximation provides high accuracy even for ti,t2 ~ 2 
[A]: 
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li A = 15.2 [eVA ], B = 24100 [eVA ] as in Ref. 
m, and ti, t2 are changing from 3.40 to 20.35 [A], then 
do(a) is in the 2.914 - 2.927[A] range. If ti and t2 tend 
to infinity then do = d[)^^ = 2.931 [A]. 

Using exact analytical formula for potential from our 
previous work2^, we can calculate the accurate value of 
equilibrium distance do{ti,t2). The maximal difference 

(2) 

between exact value do(ii,t2) and approximation do (a) 
does not exceed 0.1%. 

From Eqs. (3), (4), and (5) we have for the potential 
energy 



(p{ti,t2,d) 
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where Ci = TT^{-A/3d^ + 2B/45d^) and 
—Air'^/iSd'^ are the parameters which do not depend on 
the radii of tubes. 

Substitution of equilibrium vdW gap do into Eq. (3) 
gives us the potential well 
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FIG. 2: (a) Uniform curves for potentials and (b) forces be- 
tween two SWNTs. 



If we use the first approximation d^^ for equilibrium 

spacing (6) then Cl = -A^/^VTOtt^/QVB and = 
-744/322/351/3^2/9551/3^ 

In our previous work24 we have defined 950(^11^2) on 
the base of exact analytical formulas. After that we 
have described the potential V'o(^i!^2) by the same Eq. 
(10) and have numerically fitted = -0.19285 [eVA] 

and = -0.05847 [eVA ]. Now we have proved 
this form of dependence and have analytically evaluated 

C? = -0.1928 [eVA] and = -0.053 [eVA^], 

Comparison between exact magnitude of well depth 
and approximating value for armchair SWNTs of differ- 
ent sizes is shown in Table I. The maximal error is 3%. 

On the base of Eqs. (7), (9), and (10) it is easy to 
express the uniform curve when ti and t2 tend to infinity: 
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We note that the uniform curve depends only on di- 
mensionless distance and does not contain any material 
properties. The uniform curve ^(d) is shown in Fig. 2(a) 
as a solid blue line. The plots for SWNTs of different 
radii fall on the uniform curve with accuracy of line thick- 
ness. The splitting between all plots does not exceed 
0.013 at d = 2. 

The vdW force for two crossed SWNTs can be obtained 
by simple differentiation dip/dd of vdW potential from 
Eq. (9): 
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where Di = TT^{A/d'^ - 2B/5d^") and D2 = An'^/2Ad^. 

Using Eq. (12) we can find the distance dmax where 
the attractive vdW force reaches maximum. If ti and t2 
tend to infinity in the first approximation we have: 
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The second approximation of higher accuracy is: 
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TABLE I; Comparison of exact potential depth \(po\ [eV] (upper-right side)^ and approximation (lower-left side) for two 
SWNTs crossed at right angle 



Tube type 


(5,5) 


(10,10 


(15,15) 


(20,20) 


(25,25) 


(30,30) 


Radius [A] 


3.40 


6.79 


10.18 


13.57 


16.96 


20.35 


(5,5) 3.40 


0.761 \ 0.785 


1.060 


1.277 


1.463 


1.627 


1.777 


(10,10) 6.79 


1.039 


1.415 \ 1.433 


1.726 


1.977 


2.200 


2.403 


(15,15) 10.18 


1.256 


1.711 


2.069 \ 2.080 


2.382 


2.651 


2.895 


(20,20) 13.57 


1.442 


1.963 


2.373 


2.722 \ 2.729 


3.037 


3.317 


(25,25) 16.96 


1.606 


2.186 


2.643 


3.031 


3.376 \ 3.380 


3.691 


(30,30) 20.35 


1.755 


2.388 


2.887 


3.312 


3.688 


4.029 \ 4.031 



If we use the first approximation (13) for dma 
(12) then we have the maximal attractive force: 



in Eq. 
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where D™'^'^ 



^max 



= 3A5/3vrV5B2/3 and 
A3/27r2/24/B. 

If we define the uniform curve for vdW force as a plot 
of = -F/|iSnax| against d/ = d/dmax then we have: 
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The uniform curve for vdW force between two crossed 
SWNTs is shown in Fig. 2(b) as a solid blue line. 

Let's consider the vdW interaction between two paral- 
lel SWNTs. 

Using expansion in small parameter for integral (2) we 
get approximating formula where only the main term of 
the expansion is taken into account: 



(p[ti,t2,d) = V 
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here </> is expressed in units of energy per unit of length. 

Analogically to previous part of present work we can 
calculate the equilibrium vdW distance, the potential 
well, and the maximal attractive force for interaction be- 
tween two parallel SWNTs. 

The equilibrium distance is 
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If we use the lattice constant a = 2.49 [A], the at- 
tractive A = 15.1636 [eVA ] and repulsive constants 
B = 24052 [eVA ] as in workil then the equilibrium 
gap is do = 3.174 [A]. 

The equilibrium vdW potential per unit of length is 
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Following to workii we calculated the well depth for 
SWNTs of different radii from 2 to 22 [A]. 

Comparison between exact magnitude of well depthii 
and approximating value is shown in Table II. We can 
conclude that the maximal difference consists 6.9%. 

The equation for uniform curve is 
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The uniform curve is shown in Fig 2(a) as a dotted 
red line. The splitting between all plots for SWNTs of 
different radii does not exceed 0.015 at d = 2. Remark- 
able that the universal curve proposed by Girifalco L.A. 
et al^ 



3.41 



3.13d + 0.28 



-0.4 



3.41 



3.13d + 0.28 



is in very good agreement with our Eq. (20). If d > 1 
then ma.x\(j)Q(d) -^(d)| = 0.019 at d = 1.7. 

According to Girifalco L.A et al.— the universal curve 
is the plot of = (j)/\(j)Q\ against d = {R—ti—t2 — S)/{Ro — 
ti — t2 — S), where R is the distance between centers of 
graphitic structures, Rq is the equilibrium spacing at the 
minimum energy for the two interacting entities, ti and 
t2 are the radii of tubes or fuUerenes, and S is some pa- 
rameter. For interaction between two parallel SWNTs 
one has 6 = 0. In other cases of interaction between 
graphitic structures (Ceo — Coo, Coo-SWNT etc.) the fit- 
ting parameter S is used to adjust a plot to the universal 
curve. We prefer to not use any fitting parameter at all. 
Thus we have two different uniform curves for parallel 
and crossed SWNTs (see Fig. 2). 

The vdW force per unit of length for two parallel 
SWNTs is 



^2^2 J^ht£( 
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The distance where the attractive vdW force reaches 
maximum is 
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TABLE 11: Comparison of numerically calculated potential depth \cj>o\ [eV/A] (upper-right side)i2, and our approximation 
(lower-left side) for two parallel SWNTs 



Radius [A] 


2 


6 


10 


14 


18 


22 


2 50. 


88 \ 48.19 


61.82 


67.40 


70.53 


72.58 


74.00 


6 


62.31 8 


i8.12 \ 84.74 


95.63 


102.27 


106.76 


110.01 


10 


65.68 


98.53 


113.77 \ N/A 


119.66 


126.26 


131.17 


14 


67.31 


104.27 


122.88 


134.61 \ 131.26 


139.62 


145.97 


18 


68.26 


107.93 


129.00 


142.78 


152.64 \ 149.47 


157.07 


22 


68.89 


110.47 


133.40 


148.82 


160.09 


168.74 \ 165.76 



The maximal attractive force per unit of length is 



/maxitl ,t2)-iy 4g^8g3/4 ^3/4 \J + 
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If we define the uniform curve for vdW force between 
two parallel SWNTs as a plot of / = //l/maxl against 
df — d/dmax then we have: 
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This uniform curve is shown in Fig 2(b) as a dotted 
red line. 

In summary, we applied Lennard- Jones potential and 
method of the smeared out approximation suggested by 
L.A. Girifalco to study interaction between two SWNTs 



of different diameters. Using expansion in small param- 
eter we obtained very simple and accurate algebraic ex- 
pressions of the vdW potential and the force for two par- 
allel and crossed carbon nanotubes. It is found that inter- 
action between parallel and crossed tubes are described 
by different uniform curves. The expressions of universal 
curves contain only dimensionless distance as a parame- 
ter and do not depend on other factors. We gave the ex- 
plicit functions for equilibrium vdW distance, well depth, 
and maximal attractive force. We plotted uniform poten- 
tial curves and uniform force curves for SWNTs. These 
results may be used as a guide for analysis of experimen- 
tal data to investigate interaction between nanotubes of 
various natures. 
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